Bicuspid aortic valves undergo excessive strain during opening: A simulation study

ObjectiveThe objective of this study was to examine the influence of the morphologic characteristics of the bicuspid aortic valve on its disease progression by comparing the motion, stress/strain distribution, and blood flow of normal and stenotic tricuspid valves using simulation models.MethodsBicuspid, stenotic tricuspid with commissural fusion or thickened leaflet, and normal aortic valves were modeled with internal blood flow. Blood flow and the motion of aortic valve leaflets were studied using fluid–structure interaction finite element analysis, and stress/strain (curvature) distributions were calculated during the cardiac cycle. To mimic disease progression, we modified the local thickness of the leaflet where the bending stress was above a threshold.ResultsTransvalvular pressure gradient was greater in the bicuspid valve compared with the stenotic tricuspid valve with a similar valvular area. The bending strain (curvature) increased in both stenotic tricuspid and bicuspid valves, but a greater increase was observed in the bicuspid valve, and this was concentrated on the midline of the fused leaflets. During disease progression analysis, severity of the stenosis increased only in the bicuspid aortic valve model in terms of valvular area and pressure gradient.ConclusionsThe characteristic morphology of the bicuspid valve creates excessive bending strain on the leaflets during ventricular ejection. Such mechanical stress may be responsible for the rapid progression of this disease

Approximation for Cooperative Interactions of a Spatially-Detailed Cardiac Sarcomere Model

We developed a novel ordinary differential equation (ODE) model, which produced results that correlated well with the Monte Carlo (MC) simulation when applied to a spatially-detailed model of the cardiac sarcomere. Configuration of the novel ODE model was based on the Ising model of myofilaments, with the “co-operative activation” effect introduced to incorporate nearest-neighbor interactions. First, a set of parameters was estimated using arbitrary Ca transient data to reproduce the combinational probability for the states of three consecutive regulatory units, using single unit probabilities for central and neighboring units in the MC simulation. The parameter set thus obtained enabled the calculation of the state transition of each unit using the ODE model with reference to the neighboring states. The present ODE model not only provided good agreement with the MC simulation results but was also capable of reproducing a wide range of experimental results under both steady-state and dynamic conditions including shortening twitch. The simulation results suggested that the nearest-neighbor interaction is a reasonable approximation of the cooperativity based on end-to-end interactions. Utilizing the modified ODE model resulted in a reduction in computational costs but maintained spatial integrity and co-operative effects, making it a powerful tool in cardiac modeling